The Relation Between Cool Cluster Cores and Herschel-Detected 
Star Formation in Brightest Cluster Galaxies 



T. D. Rawle 1 , A. C. Edge 2 , E. Egami 1 , M. Rex 1 , G. P. Smith 3 , B. Altieri 4 , A. Fiedler 1 , 
C. P. Haines 1,3 , M. J. Pereira 1 , P. G. Perez-Gonzalez 5,6 , J. Portouw 1 , I. Valtchanov 4 , 
G. Walth 1 , P. P. van der Werf 7 , M. Zemcov 8,9 

trawleOas . arizona . edu 



ABSTRACT 

We present far-infrared (FIR) analysis of 68 Brightest Cluster Galaxies (BCGs) at 0.08 < 
z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry 
spanning the peak of the dust component (24-500 /im) , we calculate the obscured star formation 
rate (SFR). 22±f 2 3 % of the BCGs are detected in the far-infrared, with SFR = 1-150 M Q yr^ 1 . 
The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core 
clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is 
influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Ha 
emission is also correlated with obscured star formation. For all but the most luminous BCGs 
(Ltir > 2xlO n L Q ), only a small (<0.4 mag) reddening correction is required for SFR(Ha) to 
agree with SFRfir. The relatively low Ha extinction (dust obscuration), compared to values 
reported for the general star-forming population, lends further weight to an alternate (external) 
origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that 
the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate form 
normal stellar mass loss. 

Subject headings: Galaxies: clusters: general - Galaxies: elliptical and lenticular, cD - galaxies: star 
formation - infrared: galaxies 
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1. Introduction 

In the central regions of rich galaxy clusters, 
the intracluster medium can be sufficiently dense 
that cooling to 10 4 K occu rs on timescales shorter 
than the cluster lifetime dCowie fc Binnev 1977 : 
Fabian fc Nulsenl Il977t lEdge et al.l Il992h . How- 
ever, X-ray observations indicate that at low 
gas temperatures (a factor of a few below am- 
bient), the cooling rates are slower than expected 
for a simple run-away cooling flow model (e.g 



iTamura et al.ll200lt|Peterson fe Fabianll2006l ). Ac- 
tive Galactic Nuclei (AGN) outbursts could in- 
ject sufficient energy into the intracluster medium 
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to suppress cooling. Although the underlying 
physics is still not well understood, observation 
of radio bubbles, so called as they are regions 
relatively empty of thermal X- ray gas coinci- 
dent with the radio lobes (e.g. iBohringer et al 



1993t iMcNamara et alJkoOOl ). s upport AGN heat 



ing as a viable mechani sm (lEdge et al.l 11992 : 



McNamara fc Nulsenl 120071 ). 



Many clusters possess a large elliptical galaxy 
at the minimum of the potential well. These 
Brightest Cluster Galaxies (BCGs) are often more 
than a magnitude brighter in optical bands than 
the second ranked galaxy. In contrast to the 
majority of massive cluster ellipticals, which are 
optically red and quiescent, some BCGs contain 
significant cool gas and exhibit signs of star for- 
mation. These BCGs exhibit far- infrared atomic 
cooling lines (jMittal et al . 2011), and millime- 
ter molecular transitions (H 2 and CO emission) 
that reveal large quantities of gas at a range 
of cool temperatures (T <-, 2000 K; Edge et all 
20021 lEgami et al.l l2006bl I Johnstone et al.l 120071 



Edge et al. 2010bl ). Strong optical emission lines 



are often used as a proxy f or this molecular gas 
(e.g. Cavagnolo et al. 20081), with measured line 



ratios typical of HII regions icrawford et al"1ll999l ; 
IConselice. Gallagher fc Wvsd l200lh . Ha imaging 
of nearby BCGs appears to show emission trac- 
ing filamentary structure well beyond the central 
AG N, but apparently crea ted by the radio bubbles 
(e.g. ICrawford et "aT]l2005l) . 

The infrared dust component offers an inde- 
pendent probe of the star formation rate (SFR). 
In the Spitzer era, great sensitivity at 24 /Ltm of- 
fered the abilit y to estimate SFR fo r large samples 
of B CGs (e.g. lEgami et al] l2006at lOuillen et al 



2008 ). For instance, usin g the latter of these sam- 



ples, (0T)eaei~aD (|2008l ) found that the cluster X- 
ray emission is correlated with IR luminosity for 
Ha-detected BCGs, suggesting that cool gas orig- 
inally associated with the intracluster medium is 
fueling star formation in BCGs. This interplay is 
still not well understood, with only a small frac- 
tion of the available cold gas forming stars (< 100 
Moyr" 1 ). 

With the advent of sensitive far-infrared (FIR) 
observations from the Herschel Space Observatory, 
it is now, for the first time, possible for a direct 
measurement of the dust component in a large 
sample of BCGs. This study uses Herschel pho- 



tometry to calculate the total IR luminosity for 
BCGs in clusters at z ~ 0.08-1.0. Dust prop- 
erties and SFR can be constrained tightly, with- 
out the need to estimate from the mid-infrared. 
This is particularly important at these redshifts, 
as the varying intensity of PAH and silicate ab- 
sorption features within the 24 /j,m band causes 
non-systemmatic uncertainty in the extrapolated 
IR luminosity. 

The paper is organized as follows. Section [5] 
details the far-infrared observations, while Sec- 
tion [3] describes the important process of matching 
flux at different wavelengths (countcrparting) , and 
tabulates the final photometry. Section|4]describes 
the SED template fitting procedure and presents 
the integrated IR luminosities. In Section [5] we 
discuss the dust and star formation properties in 
the context of cool-core cluster indicators. The 
primary conclusions are summarized in Section [6] 

2. Observations 

This study employs far-infrared data from the 



ESA Herschel Space Observatory (jPilbratt et al 



2010h . Specifically, we use photometry from the 
Photod etector Array Camer a and Spectrometer 
(PACS: IPoglitsch et al.l[2010h at 70, 100, 160 fxm, 
and from the Spe ctral and Photome tric Imaging 
Receiver (SPIRE: [Griffin et al.lElloh at 250, 350, 
500 /im. These bands straddle the peak of the far- 
infrared spectral energy distribution (SED), for 
sources up to z~3. The respective beam sizes are 
larger with increasing wavelength, with FWHM = 
5.2", 7.7", 12", 18", 25", 36". 

Three Herschel key programs contribute to the 
full sample of 67 cluster fiel ds (68 BCGs): Her - 
schel Lensing Survey (HLS; lEgami et al.l l20l6T) . 
Local Cluster S ubstructure Survey (LoCuSS; 
Smith et al. l2010l ) and a study of BC Gs in known 
cool-core clusters (jEdge et al.ll2010bl E10). Her- 
schel data were processed using an augmented ver- 
sion of_the standard reduction pipelines (HIPE; 
Ott 2010f) . Major modifications are as follows. 
During PACS reduction, high-pass filtering of the 
time-stream is necessary to remove the 1/f noise 
drift. In an attempt to avoid bright source ringing, 
we implement a simple masking algorithm prior 
to filtering. In all Herschel processing, signal-to- 
noise was maximized by including data frames 
beyond the nominal science scan legs, based on 
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instantaneous scan speed: within 10% of science 
scan speed for PACS; >0.5"/s for SPIRE. In real- 
ity, turnaround frames only increase sensitivity at 
the BCG position for HLS PACS, where the total 
map size is comparable to a single array footprint. 
For further details on the Herschel data reduction, 
please refer to the individual survey papers. 

The infrared SEDs for each BCG also include, 
where available, data from the Spitzer archive 
(IRAC a nd MIPS), the 2 MAS S 2.2 /im if -band 



catalog (ISkrutskie et all l2006h and the WISE 
(jWright et all 12010) ] preliminary data release (at 
3.4, 4.6, 12, and 22 /mi). Although coverage in 
these near /mid-infrared bands is not complete, 
none of our conclusions are sensitive to their avail- 
ability. 

2.1. Herschel Lensing Survey (HLS) 

HLS ( Egami et al.ll2010l) was primarily devised 



to exploit the lensing effect of massive clusters, 
probing beyond the confusion limit of the Herschel 
instruments to obse rve intrinsically faint, high- 
redshift sources (e.g. Rex et al. 2010l ). HLS com- 
prises deep Herschel data in 5 bands, 100-500 /im. 
PACS imaging at 100 and 160 fim has mean 3d 
limits of 2.4 and 4.7 mJy, while in the three SPIRE 
bands, the typical 3 a limits are 9.4, 10.6, 12.0 mJy 
respectively. 

The survey observed 44 clusters (0.2 < z < 1.0, 
with the majority at 0.2 < z < 0.4) centered on 
the diffuse X-ray peak. More than three quarters 
of the clusters are covered by Spitzer MIPS 24 /im 
imaging, nearly all have IRAC data at 3.6 and 4.5 
/im, while approximately a third also have 5.8, and 
8 /im. Crucially for BCG identification, all of the 
cluster cores have been observed in at least one 
HST band. 

2.2. Local Cluster Substructure Survey 
(LoCuSS) 

LoCuSS is a multi-wavelength (X-ray, through 
far-UV, optical, far-IR, to radio), wide-field sur- 
vey of X-ray luminous galaxy clusters at 0.15 < z 
< 0.3. One of the main motivations is to probe 
the physics responsible for the transformation of 
galaxies from late- to early-types in the cluster in- 
fall regions and su rrounding large-scale structure 
( Smith et al.ll201ol) . A key dataset supporting this 
science goal is Herschel observations at 100-500 



/im of a representative sample of 31 clusters. The 
data covers a wider field of view (~25'x25') than 
HLS, to mean 3 a depths of 13.0, 17.0, 14.0, 18.9, 
20.4 mJy (100-500 /mi). 

The 31 cluster Herschel sample is referred to in 
this paper as LoCuSS. Eleven of the clusters are 
also in HLS, and for these, we adopt the deeper 
HLS Herschel data. The LoCuSS sample has full 
coverage in the optical from both H ST and Subaru 
Suprime-Cam (jOkabe et al j |2010l). Svitzer MIPS 
24 /zm photometry is complete ( Haines et al.l f2009. 
Haines et al. in prep.), while IRAC imaging is 
available for more than half of the clusters. MIPS 
70 /an photo metry was pub l ished for three Lo- 
CuSS BCGs in lQuillen et all (120081) . 



2.3. Known cool-core clusters (E10) 

The third component explores far-infrared line 
emission and continuum properties for a sample of 
11 well-studied BCGs in cool-core clusters. The 
BCGs were selected by optical emission and X- 
ray properties as the most likely FIR-luminous 
BCGs in the local Universe (z < 0.3). The Her- 
schel obser vations pertinent fo r this study were 
published in lEdge et al.l (|2010bl E10) and includes 
all 6 PACS+SPIRE bands (70-500 /im) for three 
BCGs (A1068, Z3146, A2597). These sources have 
extensive existing optical and infrared photom- 
etry, as well as submillimeter observations from 
JCMT/SCUBA. 

3. Multi-band photometry 

In this section, we describe the process of al- 
locating multi-band infrared flux to each optical 
BCG. This applies to the 64 clusters in the HLS 
and LoCuSS samples. Photometry for the remain- 
ing three BCGs was taken directly from E10. 

3.1. BCG identification 

For each cluster, the first task was to iden- 
tify the BCG, which we define as the optically 
brightest source in the cluster. We verified lit- 
erature positions by examining all available opti- 
cal data, including central HST imaging in every 
case. Two of the clusters are generally consid- 
ered to consist of two components, eac h with their 
own BCG: A773 dCrawford et al.lll999h and A1758 
([David &: Kempnerl 120041 ) . We confirm this selec- 



tion, and retain both sources for analysis (labelled 
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north and south). Our conclusions are indepen- 
dent of the inclusion of any combination of these 
BCGs. MACS0025.4-1222 (the "Baby Bullet") 
lacks an obvious BCG, and is removed from the 
sample. The left- most columns in Table [1] present 
the clusters and supply optical positions for the 65 
BCGs (including two pairs). 

3.2. Counterparts 

The PACS beam size (~8 arcsec) has a large 
probability of encompassing multiple optical 
sources. The most effective method of evaluating 
whether long- wavelength flux is associated with an 
optical source is to 'step through' the wavelength 
range using all available near- and mid-infrared 
data. It is important to note that a priori knowl- 
edge or expectation of whether a particular BCG 
should emit in the far-infrared was not used dur- 
ing the allocation of flux to sources. 

In all cases, the BCG is identifiable in the 
2MASS and Spitzer IRAC near-infrared bands, 
when such data is available (only 7 BCGs lack 
IRAC coverage), providing a useful link between 
optical and far-infrared wavelengths. In many 
cases, there are no Spitzer 24 fim or Herschel 
sources near the optical BCG position, and the 
BCG clearly decreases in luminosity through the 
IRAC bands (e.g. Figure ffla-fe). In a few in- 
stances, a single mid- infrared source is located in 
the immediate vicinity of the optical BCG posi- 
tion, but on closer examination, this flux obvi- 
ously originates from a secondary IRAC source 
(e.g. Figure [It). 

A few clusters contain single, point-source FIR 
counterparts to the BCG (e.g. Figure [T]i). How- 
ever, many of the Herschel-detected BCGs require 
the subtraction of secondary sources in at least the 
longer SPIRE bands where the spatial resolution 
deteriorates (e.g. Figure \Tjp-f). The techniques 
employed for extracting flux, particularly in such 
blended and confused cases, are described in the 
next section. 

3.3. Photometry 

Photometry for Spitzer and PACS, where 
blending and source confusion is minimal, was ob- 
tained using simple ape rture photometry via th e 
SExtractor package (jBertin fc Arnoutd ll996). 
Aperture corrections were applied individually in 



each band, based on instrument handbook values 
or the latest model beam profiles. 

In the SPIRE bands, background source confu- 
sion is above the instrumental noise. Photometry 
for all sources within ^3 arcmin of the BCG was 
extracted using the Iraf routine Allstar. This 
algorithm fits the known PSF to all positions in an 
input source catalog, thereby measuring photom- 
etry for all potential sources simultaneously. Any 
significant contaminating flux within a beam of 
the BCG position should be accounted for. Where 
available, MIPS 24 /im is the favored input cata- 
log, as the majority of Herschel point sources have 
visible counterparts in the MIPS image. Other- 
wise, the longest available wavelength IRAC band 
is used instead (although the correspondence be- 
tween IRAC and FIR is not as consistent). While 
residual flux from confused background sources, 
too faint for inclusion in the MIPS catalog, may 
boost the SPIRE flux, use of such an a priori cat- 
alog should minimize the effect. Statistical de- 
boosting is beyond the scope of this study. 

Table [1] presents the Herschel fluxes (100-500 
/j,m) for the HLS and LoCuSS sample BCGs. Non- 
detections are indicated by the 3 a detection limit 
of each map (note the significant difference in 
depth between HLS and LoCuSS data). Overall, 
37 of the 65 BCGs in HLS+LoCuSS are detected 
at 24 /im. A further 7 BCGs lack observations 
at this wavelength. However, in all such cases 
the slope of the IRAC photometry together with 
non-detections in PACS, suggest no significant in- 
frared component (SFR < 2 M©yr _1 ). Of the 
37 MIPS sources, 12 are detected in both PACS 
bands (there are no sources in a single Herschel 
band only). Ten of these BCGs are also detected 
by SPIRE at 250 /im, 7 at 350 fim, while two have 
significant flux at 500 /im. 

4. SED analysis and results 

In this section we analyze the full sample of 68 
BCGs (65 HLS+LoCuSS; 3 from E10), ensuring a 
homogeneous characterization of the far-infrared 
component. 

4.1. Infrared SED fitting 

Far-infrared dust properties are derived from 
the best-fitting template across all photometric 
data A bs > 24 ^m. We employ weighted fits 
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Fig. 1.— 11-band thumbnails (HST, IRAC, MIPS, PACS, SPIRE) for a representative selection of BCGs, 
identified by the cyan square (each frame is 40"x40"): (a) AS1063, (b) A1914, (c) A2744, (d) A1835, (e) 
A383, (f) A2667. See Section for details. 



to the recent templates from Rieke et al. (2009, 
R09). These templates are derived from eleven 
star-formation-dominated local (ultra-) luminous 
IR galaxies (LIRGs) with high quality data across 
the infrared wavelengths, which were combined 
and interpolated to produce a library of spectra 
for galaxies at specific luminosity. Recently, incon- 
sistencies have been discovered between the shape 
of the FIR component in these local sources and 
more distant galaxy populations. For instance, lu- 
minous sources at high redshift appear to be fit 
better by local t emplates of lower luminosity (e.g. 



Rex et alJlioioh . while a subset of cluster mem- 



ber galaxies have warmer dust than expected for 
their luminosity, an d are better fit by higher lu- 
minosity templates ( Rawle et al.l 120101 Rawle et 
al., in prep). Here, we are primarily interested in 
the total infrared luminosity, and so disregard the 
nominal luminosity class of the templates, simply 
integrating the template which best fits the ob- 
served shape of the dust component (with appro- 
priate normalization). 

Total infrared luminosity (£tir) is calculated 



by integrating the best fit template over the (rest 
frame) wavelength range Ao = 8-1000 /xm. Un- 
certainty in luminosity is estimated by repeating 
the full fitting procedure for many Monte Carlo re- 
alizations of the observed photometric data. Star 
formation rate (SFR) is derive d directly fr o m Lt ib. 
following the simple relation of lKennicutt (19980 
Figures IMTT1 in Appendix present the observed 
FIR SED for Herschel-detected BCGs, and indi- 
cates the best fit templates and subsequent values 
for Ltir and SFR. 

As a sanity check, Ltir and S FR are also de- 
rived using the template s from Chary fc Elbaj 
(|200lh and bale fc Heloul (|2002f ). The best fit- 
ting template from both are selected as described 
above for R09. SFRs derived from these tem- 
plate sets agree with the R09 values within 6% and 
2% respectively (Figure [5]) . All three template li- 
braries are unable to match the observed FIR and 
24 /jm simultaneously for a significant number of 
sources. This is most likely due to PAH and/or 
silicate absorption features within the MIPS 24 



1 SFRir [M Q yr" 1 ] = 4.5xl0" 44 L m [erg s" 1 ] 
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Fig. 2.— For FIR-bright BCGs, o bscured SFR de- 
rived using the Rieke et al. I d2(M R09) templates 
compa red to SFRs deri ved from Chary fc Elbaa 
(|200lL green circles) and bale fc Heloul (|2002i or- 
ange squares) templates, and from e xtrapolating 



the 24 /im flux via equation 14 of iRieke et al 
(blue triangles) . All ASFRs are plotted as a func- 
tion of BCG redshift, with the majority of errors 
well within the symbol size. The labeled major 
outlier (SFR 24 > SFR R09 ) is Z2089. 



/im band (rest frame Ao ~ 16-20 /im at these red- 
shifts). SFR can be estimated direct ly from the 
24 /xm flux, following equation 14 from lRieke et all 
(2009). SFR24 disagrees with the Herschel-derived 
SFR (from integrated Ltir) by an average of 20% 
(Figure [2). The major outlier (Z2089; SFR ~ 
O.O5XSFR24) is described further in subsequent 
sections. 

Characteristic dust temperature is calculated 
via a modified blackbody of the form 



S v = N^B v (T dust ) 



(1) 



where S v is flux density, j5 is dust emissivity index 
(fixed at 1.5) and B v {T dnst ) is the Planck black- 
body radiation function for a source at character- 
istic temperature Tk nst- /?Jdust remains approx- 
imate ly constant (jBlain. Barnard fc Chapman 
l2003f ). so fixing /?=2.0 instead, would systemat- 
ically lower the derived temperature (by ^3 K at 
30 K). Photometry at wavelengths shorter than 
100 /im are not used in the blackbody fit. 

Fits for two modified blackbodies, as employed 
in E10, are strongly dependent on the inclusion of 
24 /im. With the mid-infrared, best fit tempera- 
tures tend towards ^20 K and ~50 K, with much 
the same Ltir as the best fitting template. Ex- 
cluding 24 /im, the fits are less constrained and 
the two components become more similar (~30- 
40 K), or one component becomes insignificant 
and the solution tends towards the single modified 
blackbody. At most, only six of our BCG sample 
(including all three from E10) have sufficient pho- 
tometry to fit a well constrained two-component 
model. We avoid two component blackbodies in 
favor of the single characteristic temperature and 
£tir derived from the best fitting template. 

Template and blackbody fitting is only em- 
ployed for BCGs with a detection in at least two 
Herschel bands (although in practice, no sources 
have a detection in only one band). For the 
remaining BCGs, an upper limit on I/tir, and 
hence SFR, can be calculated from the Herschel 
detection limits. The detection limits are signifi- 
cantly different for each survey (but relatively sta- 
ble map-to-map within a sub-sample, e.g. 1 a de- 
viation in the HLS limits is 0.1 and 0.3 mJy for 
the two PACS bands), and the limit for Ltir is 
also a function of redshift. An additional limit 
on luminosity is estimated from the 24 /im data 
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(where available), and is used in subsequent fig- 
ures if it presents a tighter constraint than the 
Herschel-derived limit. 

4.2. Total infrared luminosities 

The upper portion of Table[2]presents Ltir and 
SFR (as derived from the R09 fits) for all Herschel- 
dctcctcd BCGs in the sample. The 24 /tin-derived 
values are shown for completeness. The majority 
of BCGs in this sample are not luminous in the 
far-infrared (Figure E). Only 15 of the 68 BCGs 
(22±|J%) are confirmed as L T ir > 10 10 L (SFR 
> 2 Moyr" 1 ), all in the redshift range 0.05 < z < 
0.6. The absolute fraction agrees approximately 
with previous studies u sing Spitzer mid-infr ared 
data alone: 2 7^% fromlEgami et all (|2006at ) and 
43±?% from lO'Dea et al.l (|2008h . The difference 
likely results from the Herschel detection limit of 
our data compared to the 2 4 /xm sensitivity, and 
fro m the Ha selec tion used bv lQuillen et al. ( 20081 ) 
and lO 'Pea et all which implicitly favors cool-core 



clusters. Seven BCGs (10±|;J%) are LIRGs (i.e. 
Ltir > 10 11 L Q ; SFR > 20 Meyr^ 1 ), which is 
again similar to, but marg inally smaller than, the 
15±5% from lO 'Pea et all The brightest BCG in 
our sample is A1835, with an Ltir = 8.5 xlO 11 
L s . 

The combined cluster sample is in no sense com- 
plete, with a bias towards the most massive and 
well studied clusters (primary targets for lensing) 
and containing three BCGs specifically selected 
due to their anticipated high FIR luminosity. The 

31 LoCuSS clusters can be considered statistically 
identical to a volume limit ed sample selected on 
X-ray luminosity only (see Smith et al.l l2010l for 
more details). The clusters within LoCuSS (in- 
cluding the 11 overlapping with HLS), contain 

32 BCGs, eight of which are detected in the far- 
infrared. This fractional value (25tg° 2 °%) is very 
similar to our full BCG sample. However, an ad- 
ditional 3 LoCuSS BCGs have a 24 /zm-estimated 
Ltir > 10 10 L Q , but fall below the Herschel de- 
tection limit. Including these, the fraction of star- 
forming BCGs in the LoCuSS sample (SFR > 2 
Moyr" 1 ) is 34™%. 

The nominal detection limits for each sample 
are also illustrated in Figure OH and they effec- 
tively vary with redshift as a power law, meaning 
the observations probe ~2.5 dex in Ltir deeper for 



the closest clusters (z^O.l) than the most distant 
(z~1.0). The Herschel observations for HLS probe 
Ltir > 10 10 L Q for all clusters z < 0.35. The 
variation of the detection limit, combined with the 
small number of detections and over-abundance of 
sources at z ~ 0.2-0.3 (due to the LoCuSS selec- 
tion criterion) makes it difficult to comment on 
possible evolution over time. 

4.3. AGN contamination 

Near- and mid-infrared studies have reported 
that IR-bright BCGs show decreasing flux through 
the IRAC band s, as expected for a stellar domi- 
nated SEP (e.g. lEgami et al.ll2006al) . In contrast, 
an AGN-dominated source would exhibit a power- 
law increase through the near-infrared. The ma- 
jority of BCGs in this present sample also appear 
to be dominated by the stellar component, and 
Ltir estimated from 24 /im does not contradict 
the Herschel detection limits. MS2137 (.2=0.31) 
is an exception to this trend, with a luminosity 
estimate from 24 /im ~0.5 dex higher than the 
nominal detection limit of the Herschel data. The 
IRAC bands suggest this is not due to a dominant 
AGN. We return to this sources in Section [5~T1 

Z2089 is the one BCG to show a clear power- 
law increase from ^5-70 /jm (in IRAC/MIPS and 
WISE photometry). The IRS spectrum (Egami 
et al., in prep.) confirms a power-law continuum 
increase for 5-35 /tm, while optical spectroscopy 
(MMT Hectospec; Pereira et al., in prep.) shows 
broad emission and line ratios suggestive of a 
strong AGN. The SEP for Z2089 shows a large dis- 
crepancy between MIPS 70 /im and the Herschel 
observations (almost a magnitude lower in Her- 
schel; see Figure HOI). IRAC, WISE and 24 fxm are 
tracking the AGN component continuum seen in 
the IRS spectrum (SFR 24 = 432 M Q yr" 1 without 
accounting for an AGN). However, it is likely that 
70 /im is substantially boosted by bright [OIII], 
[NIII] and/or [OI] emission lines (at A rost = 51.8, 
57.3, 63.3 respectively) which lie within the 70 fim 
band at z^0.24. An alternate explanation is that 
the notoriousl y shallow Spizer 70 /im imaging (in 
this case from Quillen et al. 20081) has given an er- 
roneous detection. Even if that were the case, the 
FIR peak would have to be particularly broad, or 
blue, to encompass the WISE/24 /im and PACS 
data. This highly discrepant 70 /jm point makes 
it difficult to include Spitzer data in the SEP fit- 
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as derived from the stacking analysis (see Section I5T21 for details). 
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ting procedures, and the AGN component is much 
less dominant beyond 100 /xm. Therefore, our 
Herschel-derived Ltir (^20 Mq yr _1 ) is a very 
poor fit A bs < 70 /im, but is likely to be represen- 
tative of the true obscured SFR (and is at least a 
much tighter upper limit). 

Optical emission line ratios indicate that A1835 
and to a lesser extent A1068, may have sub- 



dominant contributions fro m AGN (iQuillen et al 



20081 : iDonahue et alj|201lh . In contrast to Z2089, 
these AGN components are not identifiable from 
IRAC colors, the star-formation dominated tem- 
plates fit well at longer wavelengths, and the dust 
temperatures are reasonable given the FIR lumi- 
nosity. SFRfir is unlikely to be boosted signifi- 
cantly. 

4.4. General dust properties 

We briefly explore the general dust character- 
istics of the FIR-bright BCG population. The 
modified blackbody fit to far-infrared data re- 
veals a typical range of dust temperatures 20- 
40 K (Figure 01 • As expected for a source with 
a strong AGN component, the dust temperature 
of Z2089 is elevated relative to the normal star- 
forming population. Excluding Z2089, and ac- 
counting for the small sample size, there is no ap- 
parent evolution of BCG dust temperature with 
redshift. There is a wide range in star formation 
rates, with the three most active attaining SFR 
> 100 M yr _1 , although as previously stated, 
two of these may include a modest AGN contri- 
bution. As expected from simple models of dust 
heating, and co nfirmed by almost a decade of ob- 
servations (e.g. Blain. Barnard fc Chapmanll2003 ; 



Hwang et all2010l) . the most actively star-forming 
galaxies typically have warmer dust (grey symbols 
in Figure S]). In our BCG sample, the LIRG-type 
BCGs (Ltir > 10 11 L Q ) are on average ~5 K 
warmer than the sub-LIRGs. Generally, BCGs 
(other than the IR-bright AGN host, Z2089) are 
indistinguishable from normal star-forming galax- 
ies in this parameter space. 

5. Discussion 

5.1. Comparison to cool-core diagnostics 

Separation of cool-core and non-cool-core clus- 
ters can be difficult and often contentious. The 
X-ray cooling time is a clean primary indicator of 



cool gas in the cluster cente r dHudson et al ] l20iri 
We follow the relation from IDonahue et al.l ( 20051 ) 
for pure free-free cooling 



t c0 {K ) ~ 10 8 yrs 



K 



10keVcm z 



3/2 



kT 



X 



5keV 



(2) 

where Kq is the central entropy derived from a 
simple power law profile, and kTx is the mean 
X-ray gas tempe rature. The ACCEPT Chandra 
archival survey ([Cavagnolo et al. 2009) presents 
both Kq and kTx for 50 of the 66 clusters 
in our sample (including 12 of 14 containing 
FIR-bright BCGs). Figure [5] compares Ltir to 
this cooling time, and shows a very good anti- 
correlation for the Herschel detected BCGs (t c o < 
1 Gyr). This confirms the relation (based on 24 



Atm-ex trapolated Ltir) presented in lO'Dea et al 



(2008, igure 9), and strongly suggests that the 
obscured star formation in BCGs is connected to 
the cluster-scale cooling process. 

Ideally we would like to compare the rate 
of mass deposited (M) by the cooling central 
ICM to the condensation rate as traced by SFR. 
For a small samp le of Spitzer-detected BCGs, 
O'Dea et all (|2008h show (in their figure 10) that 



M is ^0.5-2 orders of magnitude higher than the 
SFR. This suggests that the cooling ICM is a 
plausible origin for the star formation fuel. Com- 
puting homogeneous central mass deposition rates 
for our entire cluster sample (based on e.g. AC- 
CEPT data) remains a priority, but is beyond the 
scope of this current paper. 

Section 14.31 identifies three of the most IR- 
luminous BCGs as possible hosts of mid-infrared 
bright AGN. One may worry that they are so FIR- 
bright purely because of the contribution of the 
AGN heating. However, these BCGs are also lo- 
cated in clusters with very short cooling times (< 
200 Myr), and if the correlation in Figure [5] is real, 
then they would be expected to be amongst the 
most highly star-forming. The apparent correla- 
tion (albeit very small number statistics) between 
the BCGs with the highest SFR and those with the 
stro ngest radiative AGN may be further evidence 
(e.g. lEgami et al.ll2006al) that substantial cooling 
of intracluster gas in the cluster core modifies the 
state of the AGN, part of the cyclic interplay be- 
tween cooling and heating within AGN feedback 
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(e.g. ISilk fc ReesHl998t iBower et al.ll2006h . 

We now turn to the cold molecular gas itself, 
direc tly detectable using millimeter spect roscopy 
(e.g. lEdgdliOQlt ISalome fc Combed l2003h . How- 
ever, even with the advent of Herschel spec - 
troscopy in the far- infrared (|Edge et al.l l2010al) , 
sensitivity severely restricts the sample of poten- 
tial BCG targets. In this analysis, we opt to use 
Ha line emission, a known and reliable indica- 
tor o f ~10 4 K molecular gas (e.g. ICrawford et al 
19991: ICavaenolo et al.ll2008l) . The Ha observa- 



tions are from a mix of imaging and spect r oscopi c 
observations: iDonahue. Stocke fc Gioial jl992j) 



Crawford et aD (Il999h. iQuillen et al.l (l2008h. the 
Sloan Digital Sky Survey (jAdelman-McCarthv et al 



2008) and the REFLEX BCG sample (PI: Edge) 



- + tr%) have de- 



Eighteen of the 68 BCGs (26 
tected Ha line emission (highlighted by squares in 
Figures [3] and EJ , agreeing well with the absolute 
fraction reported by previous la rge Ha surveys 



(e.g. 27±4% for 203 BCGs from ICrawford et al 
1999). The generally good correspondence be- 



tween Herschel-detection and those with mea- 
sured Ha emission (78% of Ha-detections are FIR 
sources, while 93% of FIR sources have Ha de- 
tections; Figures [3]) suggests that the two mea- 
surements are probing the same phenomenon, and 
are reaching a comparable depth. Four Ha detec- 
tions originate from BCGs undetected by Herschel 
(the properties of these are listed below the divid- 
ing line in Table H]) . Three of these sources are 
covered by the shallower LoCuSS far-infrared ob- 
servations, and Ltir(24 /im) indicates that they 
are only marginally below the Herschel detec- 
tion limit. The remaining BCG in this category 
(Z2701) has a very low L(Ha) (>0.6 dex below any 
other detection). The Ha sources undetected by 
Herschel do not, therefore, contradict a FIR-Ha 
connection. Only one BCG (A851) is detected by 
Herschel but not in Ha; we return to this source 
in the final paragraph of the section. 

There is a very good correspondence between 
Ha detection and short cooling times (Figure 
[5]). Almost all BCGs in clusters with t c o(Ko ) < 
1 Gyr are detected in Ha, confirming that the 
emission is indeed a goo d indicator of a cool- 
core cluster (as used in e.g. Cavagnolo et al. 20081: 



McDonald et al.ll201l[ ). The one notable excep- 



cooling time fits remarkably well with the 24 /im- 
estimated £tir, but the Herschel non-detection 
shows that such a luminosity is not possible (the 
limit from Herschel is ~0.5 dex lower). The BCG 
is coincident with the X-ray peak, so if it were con- 
forming to the general trend, the cold intracluster 
medium should be feeding star formation in the 
galaxy. Further explanation of this BCG requires 
additional study beyond this current paper. 

We now qualitatively compare i(Ha) and the 
total infrared luminosity (Figure [6]) . We have en- 
deavored to remove inconsistencies in L(Ha), such 
as different aperture sizes and Hq assumptions (for 
the older literature measurements) . The values are 
corrected for extinction from our Galaxy, but not 
for internal reddening, and treatment of the [Nil] 
line is different in the various surveys (see indi- 
vidual papers for more details). Despite these un- 
certainties and caveats, there appears to be a rea- 
sonable correlation (^0.34 dec scatter; grey solid 
line) over two orders of magnitude (in Ltir)- The 
four Ha-emitting BCGs undetected by Herschel 
but with Ltir(24 ^m), lie comfortably on this 
trend. On closer examination, there appears to 
be a dichotomy between the highest Ltir sources 
and the others. Fitting to only BCGs Ltir < 
2xl0 n Lq, the slope steepens significantly (black 
solid line), and the scatter decreases to 0.27 dex. 

Normal star-forming galaxies typically exhibit 
extinction A(Ha) ~ 1 mag (e.g. Kennicutt 1983t 



Hopkins et al.l 120011) . For more IR-luminous 
galaxies (Ltir > 10 11 Lq), dust may become in- 
creasingly opaque to optical wavele ngths and the 



total SFR tends towards SFRfir (iCalzetti et al 
| 20inh . From the self-consistent relations of 



iKennicuttl (1998 1 zero extinction (i.e. SFRfir 
= SFR(Ha)) would equate to log[L(Ha)] = 
log(LxiR) + 31.34 (cyan dashed line). The red 
dash-dotted line in Figure |6] indicates 1 mag of 
(uncorrected) reddening, which could easily ac- 
count for the highest luminosity BCGs. However, 
the best fit line for the lower luminosity BCGs 
(black solid) lies well above this relation. The or- 
ange dash-dot line indicates SFR(Ha) = SFRfir 
- 0.4 mag uncorrected reddening, suggesting that 
the majority of the low-luminosity sources can 
be explained by this relatively small extinction 



tion is MS2137, which is not detected by either 
Herschel or Ha, yet has t c0 (K ) ~ 140 Myr. This 



2 SFRfir [M q yr" 1 
SFR(Ho) [Mq yr" 



1.73x10- 10 L T ir [L ]; 
= 7.9xKr 42 L(Ha) [erg s^ 1 ] 
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Kennicutt Il998l) . The orange dash-dot and red dash-dotted lines show the same relation, but accounting 
for 0.4 and 1 mag respectively of uncorrected extinction in Ha flux. Corresponding SFRs, via iKennicutt 



(|1998l ) are also displayed. 
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or less. The cold gas appears to be less dusty 
than in normal star-forming galaxies, indicating a 
different origin for the fuel. 

It is worth noting that the relation indicating 
a small apparent extinction could also result from 
an intrinsically large r Ha flux, e.g. due to an ad- 
ditional heat source ([Fabian et al.ll2011l) . For this 
scenario to be plausible, the effect would need to 
be more pronounced in the less IR-bright BCGs, 
which could be true if the absolute contribution to 
the observed Ha emission was similar in all BCGs, 
and negligible in those with the highest SFRs. 

The surprising Ha non-detection of A851 sug- 
gests a much higher extinction (A(Ha) > 1 mag) 
than the typical BCG. Further examination of the 
HST optical imaging reveals a large tidal tail, 
while the BCG also has a p rojected offset from 
the X-ray peak of -280 kpc (|Bildfell et al.ll2008l ). 
10-20 x greater than the majority of this sample. 
Even if a cool-core were present it would have lit- 
tle influence. It is likely that galaxy-galaxy in- 
teraction (merger or harassment), rather than the 
cool-core phenomenon, is responsible for the star 
formation, and hence the extinction resembles the 
normal star-forming galaxy population. 

5.2. Far-infrared stacking analysis 

We now attempt to constrain the FIR luminos- 
ity of "non-cool-core" cluster BCGs. To achieve 
this, the BCGs in the HLS and LoCuSS sub- 
samples are classified into three exclusive cate- 
gories: detected in FIR and Ha, Ha-detected only 
(ie not FIR-detected), and undetected in either 
FIR or Ha. Note that (1) the three BCGs from 
E10 would all fall into the first class, and (2) we ne- 
glect a FIR-detected only category as it would con- 
sist of the debatable BCG in A851 alone. For each 
class, the Herschel maps are co-added (stacked) 
after centering on the optical BCG position and 
removing exceptionally bright nearby sources via 
point source fitting. 

Confusion noise is a time-invariant rms uncer- 
tainty and hence once reached, integrating longer 
does not improve the detection limit (in contrast 
to the more familiar instrument noise). However, 
confusion noise is still a zero-averaged gaussian 
rms and is not spatially invariant, so stacking on 
multiple known source positions from a prior cat- 
alog can detect a mean flux below the nominal 



confusion limit (e.g. Marsden et al. 20091 ). Figure 
presents the PACS and SPIRE 250 /zm stacks 
for each of three categories. For brevity, we ana- 
lyze but do not show the longer wavelength SPIRE 
stacks, which are broadly comparable to 250 /im. 

As one would expect, the FIR-detected stacks 
(11 co- added frames) exhibit a strong source in 
the center (left panels, Figure [7]). This is a useful 
verification of the stacking algorithm and optical- 
FIR alignment. The aperture corrected flux in the 
stacks equals (to within —5%) the sum of the in- 
dividual source fluxes, and the PSF is not signif- 
icantly larger than the nominal beam size for a 
single map. 

The stack of Ha-detected sources without Her- 
schel flux (center panels, Figure [7]) only comprises 
4 BCGs (3 from shallow LoCuSS maps), and does 
not significantly improve upon the HLS detection 
limit. Therefore, although 24 /jm photometry in- 
dicates that at least 3 of these sources should be 
close to the Herschel limit for an individual map, 
the stack does not reveal a detection. 

The final category, for those undetected in both 
Ha and the far-infrared, includes the vast majority 
of the sample (49 co-added frames; right panels, 
Figure [7]) . These can be considered the non-cool- 
core cluster BCGs. Again, there is no detection, 
suggesting that a BCG in this class has a mean 
3(7 upper limit in each band (100-500 /j,m) of 0.6, 
1.1, 1.3, 1.4, 1.8 mJy. These values are well be- 
low the confusio n limit of a single source in all 
band s (e.g. from iBerta et al. 2010t iNguven et aT 



2010). The dashed line in Figure [3] shows the ap- 



proximate luminosity limit (as a function of red- 
shift) derived from these Herschel flux limits. At 
any given redshift, the average FIR luminosity of 
a non-cool-core BCG is at least an order of magni- 
tude below those in cool-core clusters, suggesting 
a mean SFR < 0.17 Moyr" 1 at z - 0.2 and SFR 
< 0.42 Mgyr" 1 at z - 0.3. The stacking analy- 
sis supports a bi-modal distribution of far-infrared 
luminosity, whereby BCGs are either star-forming 
(due to a cooling ICM), or are as quiescent as nor- 
mal cluster early- types. 

5.3. Origin of the dusty gas 

We have presented strong circumstantial evi- 
dence that the fuel for star formation in FIR- 
bright BCGs originates from the cooling ICM. A 
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recent study by IVoit fc Donahue! ( 2011 ) asserts 
that the dusty gas may instead be from normal 
stellar mass loss, estimating a potential mass con- 
tribution up to 8 M Q yr _1 . We explore this further 
by calculating the stellar/dust mass ratio for our 
BCG sample, including limits on the ratio for the 
population undetected by Herschel. 

Dust masses are estimated at 500 /im via T dust 
derived from the blackbody fit, using the formula 



dust 



4^D 2 / 500 

K a bs47rB A (T dus t) 



(3) 



where /500 is the rest frame 500 /jm flux, B\ is the 
Planck function and K a bs = 0.95 cm 2 g _1 is the 



absorption coefficient at 500 /im (jDraind 120031 1. 
Rest frame 500 /im flux is estimated from the 
best fit R09 template. We derive stellar lumi- 
nosity from the 2MASS total if-band magnitude, 
adjusted for galactic extinction and K-correction. 
The redshift-dependent conversion from _R'-band 
luminosity t o stellar mass follows the formulation 
presented in lArnouts et al. l (|2007lM 

The Herschel-detected BCGs have 200 < 
M*/M dust < 6000, with a mean M*/M dust = 1570 
(Figure |8]) . The limits for individual non-detected 
LoCuSS BCGs are inconclusive as the Herschel 
imaging is fairly shallow, while those from HLS 
imaging have limits on M*/M dust of 4000-10000. 
Given that the entire stack of FIR undetected 
BCGs shows no Herschel flux (Figure [7]), we can 
bin the undetected BCGs by stellar mass (choos- 
ing 6 bins with seven BCGs each) and estimate 
the limit on M*/M dust for the stack. The limit 
for each stellar mass bin is very similar with a 
mean M*/M dust > 15000, or an order of magni- 
tude higher than the ratio for Herschel detected 
BCGs of the same stellar mass. If stellar mass loss 
were the dominant origin for the dusty gas driving 
BCG star formation, we would expect all BCGs 
of a given stellar mass to show similar star forma- 
tion rates. Otherwise, there would still be a need 
for a mechanism to trigger star formation in some 
BCGs while leaving the majority quiescent. We 
demonstrate a very good relation between SFR 
and central cooling time in Figure [5j indicating a 
much stronger probability that the star formation 
is driven and fueled by the cooling ICM. 
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Fig. 8. — Ratio of stellar mass to dust mass, 
M*/M dust , versus the stellar mass. Circular sym- 
bols represent Herschel-detected BCGs. Lower 
limits for Hershel non-detections are given by 
small colored arrows (colors indicate sample as in 
previous plots; note HLS galaxies are better con- 
strained as the limits on FIR flux, and hence dust 
mass, are lower). When FIR- undetected BCGs 
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limits in M*/M dus t can be raised to the larger 
black arrows (any stack of individually undetected 
BCGs gives a non-detection; Figure [7J- 



3 log(M/L x ) = (-0.27 ± 0.03)2 - (0.05 ± 0.03) 
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6. Summary 

We present analysis of the largest cluster sam- 
ple observed with at least 5 Herschel band (100- 
500 /urn), combining the Herschel Lensing Sur- 
vey and LoCuSS datasets. In total, the sample 
comprises 68 BCGs 0.08 < z < 1.0, of which 15 
(22^5 3%) are detected by Herschel, equivalent to 
Ltir > 10 10 L (SFR > 2 Moyr" 1 ), including 7 
classified as LIRGs. This fractional detection rate 
is similar to the LoCuSS clusters alone (25lg° 2 °%), 
which can be considered statistically identical to 
a volume-limited X-ray luminosity selected sam- 
ple. Eighteen of the 68 BCGs (26+^%) have de- 
tected Ha line emission. Generally, far-infrared 
and Ha detections corresponds well, with 24 /nm- 
based extrapolations to the far-infrared suggesting 
that most Ha-detected BCGs unseen by Herschel 
are only just below the detection limit. 

We present four associated strands of circum- 
stantial evidence which strongly suggest that the 
fuel for star formation in FIR-bright BCGs is not 
the same as in normal star-forming galaxies, but 
is directly connected to the central gas cooling in 
the cluster core. 

1 . A very good anti-correlation between Herschel- 
derived star formation rate and the X-ray 
cooling time, for BCGs in cool-core clusters 
(t c0 (K ) < 1 Gyr); 

2. FIR stacks of the Herschel-undetected BCGs 
(49 sources) show no discernible flux, plac- 
ing mean 3 a upper limits of 0.6, 1.1, 1.3, 1.4, 
1.8 mJy on the five bands (100-500 /^m), and 
implying a mean SFR for non-cool-core clus- 
ter BCGs of SFR < 1 Moyr" 1 (z < 0.4); 

3. A surprisingly low Ha extinction (~0.4 mag) 
for all but the most IR-luminous, Herschel- 
detected BCGs (L Tm < 2xlO n L Q ) com- 
pared to similar normal star formi ng-galaxies 
(^1 mag; e.g. iKennicuttl 1 19831 ). implying 
that the gas is significantly less dusty; 

4. The stellar-to-dust mass ratio of Herschel- 
detected BCGs is at least an order of magni- 
tude below the ratio for non-cool-core BCGs 
of the same stellar mass, suggesting that nor- 
mal stellar mass loss is unlikely to be the 
dominant star formation fuel source. 



We find no significant evolution in BCG dust 
properties with redshift, as previously r e porte d 



in the Ha observations of Crawford et al 



sported 

mi), 



and recently shown for all s tar-forming galaxie s 
in large survey volumes (e.g. Hwang et al. 201d() . 
However, we note that our BCG sample is not 
ideally distributed in redshift space for such an 
analysis. Despite the difference in Ha extinction, 
the FIR dust properties of BCGs explored briefly 
in this study (e.g. characteristic temperature) are 
indistinguishable from the IR-bright population as 
a whole. 

This work is partially based on observations 
made with the Herschel Space Observatory, a Eu- 
ropean Space Agency Cornerstone Mission with 
significant participation by NASA. Support for 
this work was provided by NASA through an 
award issued by JPL/Caltech. We would also like 
to thank the HSC and NHSC consortia for sup- 
port with data reduction. This publication makes 
use of data products from the Two Micron All 
Sky Survey, which is a joint project of the Univer- 
sity of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and 
Space Administration and the National Science 
Foundation. GPS acknowledges support from the 
Royal Society. 
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Table 1 

Observed Herschel far-infrared (100-500 /xm) photometry for the 65 BCGs in the HLS and 
LoCuSS samples. Non-detections denoted by the 3a limit. 



Oluster LL) 


b 

Qbcg 


Sbcg 


b 


■^cluster 


U/T C 

ri / Li 


/ioo 


/l60 


/250 


7350 


7500 




















(mjyj 


(mjyj 


(mjyj 


(mjyj 


(mjy) 


A2744 


00 


14 


20.7 


-30 


24 


00 


0.308 


H 


<2.4 


<4.5 


<7.8 


<9.2 


<10.7 


MACS0018 


00 


18 


33.6 


+16 


26 


15 


0.541 


H 


<2.3 


<4.8 


<10.0 


<10.4 


<11.4 


A0068 


00 


37 


06.9 


+09 


09 


24 


0.255 


HL 


<2.5 


<5.0 


<12.9 


<13.0 


<13.3 


A2813 


00 


43 


24.7 


-20 


37 


41 


0.292 


H 


<2.4 


<4.7 


<8.3 


<10.0 


<11.9 


A0115 


00 


55 


50.6 


+26 


24 


37 


0.197 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


Z0348 


01 


06 


49.4 


+01 


03 


22 


0.254 


L 


81.1 ± 2.8 


54.7 ± 5.8 


35.0 ± 5.2 


12.0 ± 4.3 


<20.4 


A2895 


01 


18 


11.1 


-26 


58 


12 


0.227 


H 


<2.4 


<4.7 


<7.7 


<9.1 


<11.1 


A0209 


01 


31 


52.5 


-13 


36 


41 


0.206 


HL 


<2.4 


<4.7 


<7.8 


<9.5 


<11.0 


R0142 


01 


42 


03.4 


+21 


31 


15 


0.280 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


A0267 


01 


52 


41.9 


+01 


00 


27 


0.231 


HL 


<2.5 


<4.8 


<8.0 


<9.6 


<11.2 


A0291 


02 


01 


43.1 


-02 


11 


47 


0.196 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


RCS0224 


02 


24 


34.2 


-00 


02 


32 


0.773 


H 


<2.3 


<4.7 


<8.4 


<9.6 


<10.9 


A0368 


02 


37 


27.8 


-26 


30 


29 


0.220 


H 


<2.5 


<4.6 


<7.7 


<8.9 


<10.2 


A0383 


02 


48 


03.4 


-03 


31 


44 


0.187 


HL 


12.3 ± 0.1 


32.5 ± 0.5 


14.9 ± 6.3 


14.6 ± 6.0 


<10.6 


A3084 


03 


04 


03.9 


-36 


56 


27 


0.219 


H 


<2.4 


<4.4 


<7.7 


<9.1 


<10.7 


A3088 


03 


07 


02.1 


-28 


39 


58 


0.253 


H 


6.9 ± 0.2 


6.1 ± 0.5 


<7.0 


<8.5 


<10.1 


MACS0451 


04 


51 


54.6 


+00 


06 


18 


0.430 


H 


<2.3 


<4.2 


<13.0 


<15.6 


<19.5 


A0521 


04 


54 


06.9 


-10 


13 


24 


0.253 


H 


<2.3 


<4.5 


<9.4 


<10.2 


<11.9 


AS0592 


06 


38 


45.2 


-53 


58 


22 


0.222 


H 


<2.4 


<4.7 


<15.3 


<15.7 


<16.3 


MACS0647 


06 


47 


50.7 


+70 


14 


54 


0.591 


H 


<2.3 


<4.8 


<14.7 


<14.2 


<14.5 


BULLET 


06 


58 


38.1 


-55 


57 


26 


0.296 


H 


<3.1 


<5.8 


<8.2 


<9.6 


<12.4 


MACS0717 


07 


17 


37.2 


+37 


44 


23 


0.546 


H 


<2.3 


<4.3 


<10.5 


<11.1 


<11.7 


A0586 


07 


32 


20.3 


+31 


38 


01 


0.171 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


Z1432 


07 


43 


23.1 


+17 


33 


42 


0.111 


L 


<13.0 


<17.0 








MACS0744 


07 


44 


52.8 


+39 


27 


25 


0.698 


H 


<2.3 


<4.6 


<7.9 


<9.5 


<10.7 


A0611 


08 


00 


56.8 


+36 


03 


24 


0.288 


H 


<2.3 


<4.5 


<8.8 


<9.6 


<10.7 


Z1693 


08 


25 


57.8 


+04 


14 


48 


0.225 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


A0665 


08 


30 


57.3 


+65 


50 


32 


0.182 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


A0689 


08 


37 


24.6 


+14 


58 


22 


0.279 


L 






<14.0 


<18.9 


<20.4 


Z1883 


08 


42 


55.9 


+29 


27 


27 


0.194 


L 


15.3 ± 2.2 


15.4 ± 2.9 


10.2 ± 5.2 


<18.9 


<20.4 


A0697 


08 


42 


57.6 


+36 


22 


00 


0.282 


HL 


<2.3 


<4.2 


<8.2 


<10.0 


<11.7 


Z2089 


09 


00 


36.9 


+20 


53 


40 


0.235 


L 


55.5 ± 2.5 


24.3 ± 2.9 


12.4 ± 5.2 


5.8 ± 4.3 


<20.4 


A0773S 


09 


17 


53.4 


+51 


43 


39 


0.217 


H 


<2.4 


<4.3 


<8.0 


<9.6 


<11.2 


A0773N 


09 


17 


53.5 


+51 


44 


01 


0.217 


H 


<2.4 


<4.3 


<8.0 


<9.6 


<11.2 


A0851 


09 


42 


57.5 


+46 


58 


50 


0.407 


H 


23.7 ± 0.3 


25.3 ± 0.7 


20.5 ± 5.5 


12.3 ± 6.2 


<11.0 


A0868 


09 


45 


26.4 


-08 


39 


07 


0.153 


H 


<2.4 


<5.2 


<7.7 


<8.9 


<10.6 


Z2701 


09 


52 


49.2 


+51 


53 


06 


0.214 


H 


<2.4 


<4.7 


<7.1 


<8.2 


<9.8 


A0963 


10 


17 


03.6 


+39 


02 


49 


0.206 


HL 


<2.5 


<4.9 


<7.9 


<10.0 


<11.2 


MACS1149 


11 


49 


35.7 


+22 


23 


55 


0.544 


H 


<2.4 


<4.4 


<7.6 


<9.0 


<10.5 


A1413 


11 


55 


17.8 


+23 


24 


16 


0.143 


H 


<2.4 


<4.6 


<7.7 


<8.9 


<10.6 


A1689 


13 


11 


29.5 


-01 


20 


28 


0.183 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


A 1703 


13 


15 


05.3 


+51 


49 


02 


0.281 


H 


<2.4 


<4.5 


<7.9 


<9.4 


<11.8 


A1758N 


13 


32 


38.4 


+50 


33 


35 


0.279 


HL 


<2.4 


<4.9 


<8.0 


<9.4 


<11.3 


A1758S 


13 


32 


52.1 


+50 


31 


34 


0.279 


HL 


<2.4 


<4.9 


<8.0 


<9.4 


<11.3 
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Table 1 — Continued 





b 

CtBCG 


x b 

OBCG 


^cluster 


H/L c 


J 100 


J 160 


f 

J250 


f 

J350 


f 

J 500 












(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mJy) 


A 1763 


13 35 20.1 


+41 00 04 


0.228 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


A1835 


14 01 02.1 


+02 52 43 


0.253 


L 


315.1 ± 1.8 


322.6 ± 2.9 


122.0 ± 5.2 


45.9 ± 4.3 


21 5 + 6 2 


MACS1423 


14 23 47.9 


+24 04 42 


0.543 


H 


13.0 ± 0.1 


27.7 ± 0.3 


20.8 ± 4.6 


14.4 ± 5.8 


<10.2 


A1914 


14 25 56.6 


+37 48 59 


0.171 


HL 


<2.5 


<5.0 


<7.9 


<9.0 


<11.0 


Z7160 


14 57 15.1 


+22 20 35 


0.258 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


A2218 


16 35 49.1 


+66 12 44 


0.171 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


A2219 


16 40 19.8 


+46 42 42 


0.228 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


RXJ1720 


17 20 10.1 


+26 37 32 


0.164 


HL 


9.3 ± 0.1 


7.2 ± 0.9 


<13.6 


<14.5 


<16.2 


A2345 


21 27 13.7 


-12 09 46 


0.176 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


MACS2129 


21 29 26.2 


-07 41 29 


0.589 


H 


<2.4 


<5.3 


<29.7 


<27.7 


<27.2 


RXJ2129 


21 29 40.0 


+00 05 21 


0.235 


HL 


3.9 ± 0.1 


6.5 ± 0.3 


7.1 ± 1.8 


<10.1 


<12.0 


MS2137 


21 40 15.1 


-23 39 40 


0.313 


H 


<2.4 


<4.6 


<8.0 


<9.3 


<10.3 


A2390 


21 53 36.8 


+17 41 45 


0.233 


L 


17.5 ± 2.4 


38.6 ± 5.7 


31.4 ± 5.2 


20.4 ± 4.3 


8.2 ± 6.2 


A2485 


22 48 31.1 


-16 06 26 


0.247 


L 


<13.0 


<17.0 


<14.0 


<18.9 


<20.4 


AS1063 


22 48 44.0 


-44 31 51 


0.348 


H 


<2.3 


<4.3 


<8.0 


<9.9 


<11.5 


AS1077 


22 58 48.3 


-34 48 07 


0.312 


H 


<2.4 


<4.6 


<8.0 


<9.2 


<10.0 


A2537 


23 08 22.3 


-02 11 33 


0.295 


H 


<2.5 


<5.1 


<11.9 


<12.8 


<12.6 


RCS2319 


23 19 53.4 


+00 38 13 


0.897 


H 


<2.3 


<4.5 


<7.7 


<9.0 


<10.4 


RCS2327 


23 27 27.6 


-02 04 37 


0.700 


H 


<2.5 


<4.8 


<9.6 


<11.0 


<13.5 


A2631 


23 37 39.7 


+00 16 17 


0.273 


H 


<2.4 


<4.7 


<8.2 


<9.5 


<10.8 


A2667 


23 51 39.4 


-26 05 03 


0.230 


H 


20.4 ± 0.1 


27.4 ± 0.4 


10.6 ± 2.6 


<9.5 


<11.9 


a Cluster name, including A/B designat 


ion (as in 


the literature) for multiple candidate BCGs in a 


single cluster 




b Position of the BCG in 


optical (HST) imaging 














c Subsample: 


H=HLS; L= 


=LoCuSS. For clusters in 


t both samples, HL is 


shown for reference, although the deeper HLS Herschel 



images have provided the photometry. 
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Table 2 

Derived SED parameters for Herschel and/or Ha detected BCGs. 



Cluster ID a 


H/L/E b 


£tir(24) c 


SFR24 C 




SFR F i R d 


L(Ra) e 






(L©) 


(M yr" 1 ) 


(Lo) 


(M yr- 1 ) 


log (erg s" 1 ) 


7flO A O 

Z0348 


T 

L 


2.9 ± U.U f X 11) 


50.5 ± 1.1 


1 n _i_ n nc w 1 nil 
1.9 ± 0.06X 10 


32.6 ± 1.2 


a o c n 
42.50 


A0383 


HL 


9.4 ± 0.51 xlO 9 


1.6 ± 0.1 


2.3 ± O.lOxlO 10 


4.0 ± 0.2 


41.85 


A3088 


H 


1.3 ± 0.03 xlO 10 


2.2 ± 0.1 


1.8 ± 0.06xl0 10 


3.0 ± 0.1 


41.30 


Z1883 


L 


<1.3xl0 10 


<2 


2.4 ± 0.32 xlO 10 


4.1 ± 0.5 


42.14 


Z2089 


L 


2.0 ± 0.12 xlO 12 


346.4 ± 20.6 


1.2 ± 0.06 xlO 11 


20.2 ± 1.0 


42.50 


A0851 


H 


5.7 ± 0.19 xlO 10 


9.9 ± 0.3 


1.9 ± 0.09xlO n 


32.7 ± 1.6 




Z3146 


E 


3.5 ± 0.34 xlO 11 


61.0 ± 6.0 


5.4 ± 0.29 xlO 11 


93.1 ± 5.1 


42.55 


A1068 


E 


1.5 ± 0.04 xlO 12 


254.4 ± 6.8 


5.8 ± 0.03 xlO 11 


99.3 ± 0.5 


41.94 


A1835 


L 


1.2 ± 0.01 xlO 12 


202.0 ± 0.5 


8.5 ± 0.05 xlO 11 


146.5 ± 0.8 


42.20 


MACS1423 


H 


1.4 ± 0.04 xlO 11 


24.8 ± 0.7 


2.7 ± 0.17 xlO 11 


46.0 ± 2.9 


42.00 


RXJ1720 


HL 


1.2 ± 0.19 xlO 10 


2.0 ± 0.3 


9.3 ± 0.27 xlO 9 


1.6 ± 0.1 


40.88 


RXJ2129 


HL 


6.0 ± 0.25 xlO 9 


1.0 ± 0.1 


1.0 ± 0.05 xlO 10 


1.7 ± 0.1 


40.75 


A2390 


L 


5.1 ± 0.03 xlO 10 


8.8 ± 0.1 


4.9 ± 0.49 xlO 10 


8.5 ± 0.7 


41.69 


A2597 


E 


1.2 ± 0.12 xlO 10 


2.1 ± 0.2 


1.8 ± 0.05 xlO 10 


3.1 ± 0.1 


41.48 


A2667 


H 


1.8 ± 0.02 xlO 11 


31.2 ± 0.3 


5.0 ± 0.14 xlO 10 


8.7 ± 0.2 


41.60 


A0115 


L 


1.8 ± 0.04 xlO iU 


3.0 ± 0.1 


<2.5xl0 lu 


<4 


41.03 


A0291 


L 


1.8 ± 0.04 xlO 10 


3.2 ± 0.1 


<2.5xl0 10 


<4 


41.36 


Z2701 


H 


<4.2xl0 9 


<1 


<6.4xl0 9 


<1 


40.64 


Z7160 


L 


2.4 ± 0.04 xlO 10 


4.2 ± 0.1 


<4.2xl0 10 


<9.6 


41.40 



BCGs listed beneath the dividing line were not detected in any of the FIR Herschel bands, but do have a 

measured Ha emission line flux 

a Cluster ID as in Tabled] 

b Subsample: H=HLS; L=LoCuSS; E= lEdge et all (|2010bl ). Clusters in both HLS and LoCuSS are indi- 



cated, although the deeper HLS Herschel provided the photometry. 
c Estimated from the 24 /xm photometry (via Rieke et al. 2009) 



d Derived from the best fit I Rieke et al. (2009) template 
Titerature Ha luminosity 
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A. Infrared SED fits 



Figures IMTTI present the infrared (observed frame A = 2-2000 /iin) spectral energy distribution for the 15 
Herschel-detected BCGs in HLS, LoCuSS and E10. 
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Fig. 9. — Far-infrared SEDs for Herschel-detected BCGs (sample, name and redshift shown in black). PACS 
and SPIRE detections displayed as orange and red circles respectively. Upper limits are shown as 'v' at 
the 3 cr detection limit. Additional photometry is also displayed: 2MASS A-band (magenta); Spitzer IRAC 
(grey) and MIPS (blue); WISE (green); IRAS (cyan); SCUBA (black). A modified-blackbody is fit to the 
Herschcl points (dashed line), with the temperature and estimated dust mass shown in the low er-right corner. 
Best fi t templates f rom three librar i es are also plotted: orange = iRieke et al. (2009;), cyan = Dale fc Helou 
(2002j) and green = ICharv fc Elbad ( 20011 ). Derived Z/tir and SFRfir for each are given in matching colors. 
MIPS 24 /im-dcrivcd SFR is shown in blue in the upper-right corner. 
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Fig. 10. — Far-infrared SEDs for Herschel detected BCGs. Layout as in Figure [9] 
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Fig. 11. — Far-infrared SEDs for Herschel detected BCGs. Layout as in Figure [9j 
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